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Abstract 



This paper is devoted to the spectral theory of the Schrodinger operator on 
the simplest fractal: Dyson's hierarchical lattice. An explicit description of the 
spectrum, eigenfunctions, resolvent and parabolic kernel are provided for the un- 
perturbed operator, i.e., for the Dyson hierarchical Laplacian. Positive spectrum 
is studied for the perturbations of the hierarchical Laplacian. Since the spectral 
dimension of the operator under consideration can be an arbitrary positive num- 
ber, the model allows a continuous phase transition from recurrent to transient 
underlying Markov process. This transition is also studied in the paper. 
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(N 1 Introduction 

The spectral theory of the fractals, which are similar to the infinite Sierpinski gasket (i.e. 
the spectral theory of the corresponding Laplacians) is well understood (see [6], [7], [5]). 
It has several important features: the existence of a large number of eigenvalues of infinite 
multiplicity, pure point structure of the integrated density of states, compactly supported 
eigenfunctions. These features manifest themselves in the unusual asymptotics of the heat 
kernel, the specific structure of the corresponding ^-function, etc., see [8]. 

The next natural step in the spectral theory is to study Schrodinger type operators, 
i.e., fractal Laplacian perturbed by a potential. There are two possible directions for such 
a development: analysis of the random Anderson Hamiltonians (the potential is stationary 
in space) or the study of the classical problem on the negative spectrum when the potential 
vanishes at infinity. For the first direction, see [9], [2], [10]. We will concentrate on the 
second problem in a particular case of the simplest fractal object: Dyson's hierarchical 
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Laplacian perturbed by a decaying potential. Our goal is to prove the Cwikel-Lieb- 
Rozenblum (CLR) estimates for the number of negative eigenvalues and estimates for 
Lieb-Thirring (LT) sums. These estimates depend on the spectral dimension Sh of the 
fractal (which can take an arbitrary positive value.) The most important part of the 
paper is the analysis of the spectral bifurcation near the critical dimension Sh = 2. 

The authors are very grateful to E. Akkermans, J. Avron and A. Teplyaev for useful 
discussions. 

2 Hierarchical lattice and Laplacian 

The concept of the hierarchical structure was proposed by F. Dyson [1] in his theory of 
1-D ferromagnetic phase transitions. There are several modifications of the hierarchical 
Laplacian (see [2]). We will study the simplest one, which is characterized by an integer- 
valued parameter v > 2 and a probabilistic parameter p e (0, 1). 

Description of the model. Consider a countable set X and a family of partitions 
IIo C III C II2 C ... of X (we write n r C II r .+i to mean that every element of n r is a 
subset of some element of II r+ i). The elements of n are the singleton subsets of X. They 
are denoted by Q^p and called cubes of rank zero. Each element of IT (cube of rank 

one) is a union of v different cubes of rank zero, i.e., X = UQ[ , \Qi \ = v (see Figure 
(2) 

1). Each element Q\ of 1I2 (cube of rank two) is a union of v different cubes of rank 
one, i.e., X = UQ[ 2 \ \Q\ \ = is 2 , and so on. The parameter v > 2 is one of the two basic 
parameters of the model. 




Q« Q< 2) Q< 2> 

Figure 1: An example of a hierarchical lattice with X = Z and v = 2. 

Each point x belongs to an increasing sequence of cubes of each rank r > which we 
denote by Q {r) (x), i.e., x = Q (0) (x) C Q {1) (x) c Q {2) (x) c • • • . 
The hierarchical distance dh(x,y) on X is defined as follows: 

d h (x,y) = min{r : 3Q ( f ) 3 x,y}. (1) 

We assume the following connectivity condition holds: for each x,y G X, the cubes 
Q( n \x) contain y when n is large enough, i.e., dh(x,y) < 00. 
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Note that for arbitrary z G X, dh(x,y) < max {dh,(x, z), dh(y, z)}, i.e., dh(-, •) is a 
super-metric which implies that 

p(x,y) = pf,(x,y) = ePU x ri-l, > 0, 

is also a metric. We will use it in the form 

p(x,y)=[—) -1, (2) 



i.e., (3 = In -i=. Here p6 (0, 1) is the second parameter of the "Laplacian" A^ (see formula 
(ED below). 

Now we denote by I 2 (X) the standard Hilbert space of square summable functions 
on the set X and define a self-adjoint bounded operator (the hierarchical Laplacian) 
depending on the parameter p G (0, 1) : 



A h i/> (x) = J2 



oo 
r=l 



ip [x) 



where a r = (1 — p)p r 1 , ^^ a r = 1- 



r=l 



(3) 

The random walk on (X, dh) related to the hierarchical Laplacian has a simple struc- 
ture. It spends an exponentially distributed time r (with parameter one) at each site 
x. At the moment r + it randomly selects the rank A; of a cube k > 1, with 

P{A; = r} = a r and jumps inside of Q^ k \x) with the new position x' G Q^ fc ^(x) being 
uniformly distributed. 

It is clear that A^ = A* h , Ah < 0, Sp(A/J G [—1, 0]. The following decomposition will 
play an essential role. Denote by Ik{x) the indicator function of a set K G X, i.e., Ik — 1 
on K, Ik = outside of i^. Then, for each y G X, 



OO ✓ 

W = E 

fe=i ^ 



pk—l yk 



(4) 



The validity of OH) is obvious. It is important that each term on the right is an eigen- 
function of Ah and the k th term belongs to the eigenspace Lj, defined in the following 
proposition. 

Proposition 2.1. (a) The spectrum of Ah consists of isolated eigenvalues = —p k ^ 1 , k = 
1,2,..., each of infinite multiplicity, and their limiting point A = 0. 

(b) The corresponding eigenspaces C 1 2 {X) have the following structure: Fork = 1, 



L 1 = {^G/ 2 (X): V ^(x) = for each Q^eU,} 



For k > 1, the space consists of all ip G I (X) which are constant on each cube Q\ 



>(*) 

'i 

(c) The following decomposition holds: l 2 (X) = ©^Z/r 



and have the property that YI^qW ip(x) = for each Q\ G IT fc . 
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Indeed, one can easily check that the space Lk, defined above, consists of eigenfunctions 
with the eigenvalue A& = — p k ~ 1 , and for each y G X, the k th term in (jlj) belongs to 
Thus (jlj) immediately implies (c) which justifies (a). 

Let us note that each eigenspace Lk has an orthogonal basis of compactly supported 
eigenfunctions. Such a basis in L\ consists of functions which are zero outside of a fixed 
cube and such that Yu x eQ^ ^{ x ) = 0- There are v — 1 orthogonal functions with 

the latter property for each cube Q\ . The orthogonal complement of L\ consists of 
the functions ip G ^ 2 (^0 which are constant on each cube of rank one. The basis in L 2 
is formed by functions supported by individual cubes of rank two such that ip(x) = 
on sub-cubes Q^p of rank one, and ^ = 0. One needs to specify q to guarantee the 
orthogonality of the elements of the basis. The basis in k > 1, is formed by functions 
which are supported by individual cubes of rank k and which are constant on sub-cubes 
of rank k — 1 with the sum of those constants being zero. 

Let's find the density of states for A^ and the spectral dimension Sh- We fix xq G X 
(the origin) and a positive integer N . Consider the spectral problem 

-A^ = A^; V = 0onX\gW(x ). 

(now it is more convenient to work with — A^ instead of A/J. It is easy to see (compare 
to Proposition 12. ip that the problem has the following eigenvalues: 

Ao,jv = 1 with multiplicity v N ~ x {y — 1) 
Ai = P with multiplicity v N ~ 2 {y — 1) 



Ajv-i,iv = p^ 1 with multiplicity (i/ — ly 
A at at = with multiplicity 1. 



This implies the following relation for 



A^v(A) = ^#{A M < A}. 



Proposition 2.2. As N — > oo, 



k>0:p k <\ 

where fco(A) = min{/i; > : p k < A}. Furthermore, 



dNjX) ( 1 



5l(A) + MA) + V(A) + . 
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Proposition 2.3. As A 4- 0, 



N (A) x A Sh/2 , 



ln(l/p) 

or, more precisely, 

N (A) ~ A Sh/2 /i 



In A 

hip 

/or a positive, periodic function h(z) = u^ 1 ^^ = h(z + 1). Here, {z} is the fractional 
part of a number zGK. 

The latter proposition is a consequence of the following simple calculation. If [z] is 
the integer part of z G R, then 



In A 

hip 



We will call the constant Sh = j^fj^ the spectral dimension of the triple (X, d n (-, ■), Ah). 
Transition probabilities and the resolvent for A h . Let p (t, x, y) = P x {x (t) = y} 
be the transition function of the hierarchical random walk x (t), i.e., 

dp 
dt 

and let 



Ap, p(0,x,y) = S y (x), 



OO 

-xt 



R\(x,y)= e xt p(t,x,y)dt, A > 0. 



The functions p and R\ define the bounded integral operators 

(p t f){x) = 5>(t,z,y)/(i/), 

y ex 

(R x f)(x) = J2 R ^ x >y)f(y) 

yex 

acting in Z°° (X) and I 2 (X), respectively. 

Formula (j4]) (where each term on the right is an eigenfunction of A/J and the Fourier 
method lead to the following statement 

Proposition 2.4. The transition kernel p (t , x , y) has the form: 



p{t,x,x) = ( 1 - 1 



v 



e"* + + ••■ + — ^ + - 

v v K 



e - P t e -p k t 



for each x G X. 



e' pr ~ H ( 1\ fe~ prt e~ pr+H 
Here, r = d h (x, y) is the minimal rank of the cube (x), containing the point y (see([T]) 
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Similar formulas for R\(x, y) can be obtained from (J3J) or (easier) from the proposition 
above (by integration in t): 

Proposition 2.5. For any Sh > 0, A > 0, 

R x (x , x) = - (A - pr _^ r + (l - -J ^ x + pr y. + (A + pr+ i )z/r+ i +•■ 
when r = dh(xo,x) > 0. If xq = x, then (independent of x G X), 



R\(x,x) =11 



1 1 1 



1/J[X + 1 (X+p)u (\+p s )v s 

Corollary 2.6. (a) If pu > 1 (s A = fjf- > 2), tfien /or eac/i iGl 



(6) 



Ro(x, x) = I p(t, x, x)dt = ( 1 - - J ( 1 + — + — ^ H ] = ^ — ^ < oo. 

J \ vj \ pv (puy J pi; -I 

If pv < 1 (i.e., Sh = in(i/p) — 2); ^en lim,v->+o R\(x, x) = oo. TTras £/ie random walk x(t) 
with the generator Ah is transient for Sh > 2 and recurrent for Sh < 2. 
(6) If Sh > 2 and p(xo,x) —¥ oo (see (TJj],), i/ien 

^0(^0^) = ( — - ^T^j + (pr+l^r+l - p^T^) + " ' 

1 — p c pz/(l — p) 

c - 



(pv) r " 1 (pv — 1) p s,i_2 (xo, x) ' pi/ — 1 

27ms zs one more indication of a similarity between Ah and the lattice "L d Laplacian. 

Now let's find the asymptotics of p(t, x, x) as t — > 00. The asymptotics will play an 
essential role in the spectral theory of the Schrodinger operator H = —Ah + V(x). 

Proposition 2.7. For arbitrary spectral dimension Sh, 

1 

and there exists a positive periodic function h\(z) = hi(z + 1) such that 



hi 



lnt 



p(t,x,x) = X ^J 2 P J (1 + o(l)) ast^oo. (7) 
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— s t 

Proof. The index of the maximal term in the series pit, x, x) — (1 — -) Yl^Lo e has 
order s = 0( ln |"^ ) when t — v oo. We put = [ ln nL) ] and change the order of terms in 
the series representation of p, first taking the sum over s > k and then taking the sum 
over s < k: 



p{t,x,x) =(!--) [ — + + + + 







( e~P kt 










e -P k t ' 
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- 





e p k t(l-p) e p k t(l~p 2 ) e p k t(l-^) e^ 1 -^) 

1 + + 5 + ••■ + ^ + 5 — + •• 



The relation ^^jf-^y = k + {n^y} implies that 

i f In t i 1 In t i r In t i \ f lnt i 

xT£ = r> tln (Vp) and — = e mvp) // ^RiM' = v ^RiM' 

^ ^ z/fc phi" 1 

We substitute the latter relations into ()8]) and note that {x} is a periodic function of x 
with period one. 

This and would lead to ([7j) with zero reminder term if both series in square brackets 
in ([H]) had infinitely many terms. Since the second part in the square brackets has only k 
terms we obtain (J7J) with an exponentially small reminder. □ 

The next statement provides the asymptotic expansion of R\(x, x) as A — > +0. We 
restrict ourselves to the more difficult and important case where Sh < 2. As in the previous 
proposition, the main term of the expansion contains a periodic function. We will use an 
alternative approach to show that. 

Proposition 2.8. If Sh < 2, then 

R x (x,x) = \- a u(^)+c + O(\), A^+0, « = 1 -i^ = 1 -f> 

where Co = p ^~^ is a constant and u{z + 1) = u{z) is a positive periodic function with 
period one. 

Proof. From series representation ([6]) it follows that 

R P \ R\ 



pv v(p\ + 1) 
We put R x = c + /(A). Then 

/(PA) - 1/(A) = i^A. 
pv v\p\ + 1) 
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After the substitution /(A) = A a g(X) we arrive at 



g(p\)-g(\) = C{\), C(A) 



P 2 (l-^) A l + a 



p\ + 1 



(9) 



The estimate |C(A)| < C|A 1+a |, A > 0, is valid for the function ( (this estimate was 
the goal of the subtraction of the constant cq from R\ made above). Hence the series 
fi'par = X]o°C(p^)) A > 0, converges, has order 0(A 1+a ) as A — > +0 and is a partial 
solution of equation Any solution of the homogeneous equation OH]) is a periodic 



function of ln p A 



In A 

lnp 



with period one. This completes the proof. 



□ 



Remark. The statement of the proposition and its proof remain valid if A — > in the 
complex plane, and | argA| < 37r/4. 

We conclude this section by defining two functions, 9(t) and ^(z), which are the ana- 
logues of the corresponding classical 1-D functions: 



0(t) 



-At 



dN(X) 



-pt 



-p 2 t 



+ 



+ 



+ 



r z) 



t z - l 6{t)dt 



1 

1 - - 

v 



oo 



r=0 



1 

1 - - 

V 



p V 



p z v — 1 



The formula for ^ (z) is obtained for Re z G (0, 5) with a small enough 5 > (p Kcz v > 1) 
and understood in the sense of the analytic continuation for other z. The function q has 
no complex zeros, but (compare to (8j) has infinitely many poles at z = z n 



Sh_ _|_ inn 
2 lnl/p' 



3 Elements of the general lattice spectral theory 

The functions p(t,x,y) and R\(x,y) play a central role in the analysis of the positive 
spectrum of the hierarchical Schrodinger operator 

H = A h + V(x); V>0. (10) 

With only weak assumptions on V, the positive spectrum A n = \ n (H) > of if is 
discrete (possibly, with accumulation at A = 0). Our goals are to find upper bounds on 
N (V) = #{A n > 0} and on the Lieb-Thirring sums ^ 7 (V") = E„0n) 7 ,7 > 0. Below, 
we will provide several estimates on iVo and S 1 which are valid [3],jl] for general discrete 
operators and for the operator (fTUl) in particular. 

Let X be an arbitrary countable set and let H be a bounded self-adjoint operator on 
l 2 (X) given by 

y-.y^x 
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h(x,y) = h(y,x) > for x ^ y, h(x,y) < Co < oo. 

y.y^x 

It is clear that H = H*, H < 0, \\H \\ < 2C . 

Let p(t,x,y) = P x (x(t) = y) be the transition kernel of the continuous time Markov 
chain x(t) generated by H . Of course, 



dp 
dt 



— = H p, p(0, x, y) = 8 y {x). 



We assume that x(t) is connected which means, since its time is continuous, that 
p(t, x,y) > for arbitrary x, y G X and t > 0. 

The bounds for the eigenvalues of H depend essentially on whether the process x(t) 
is transient or recurrent. If J °° p(t,x,x)dt < oo for every x G X, then x(t) is transient, 
i.e., P-a.s., x(t) — > oo as t — > oo. If J °° p(t,x,x)dt = oo for every x G X, then x{t) 
visits each state x G X infinitely many times P-a.s. and the process is called recurrent. 
It is a well-known fact that, if the chain is connected, the convergence or divergence of 
Jo°° P^Pi x i y)dt ls independent of x, y. 

Theorem 3.1. (General CLR estimate for discrete operators) . If L p(t, x, x)dt < oo, 
then for any a, a > and some c\(a), 



N (V) < #{x G X : V(x) >a} + Cl (a) ^ V(x) J p(t 



x:V(x)<a 

V{x) 



Theorem 3.2. (LT estimate). If J °° p(t,x,x)dt < oo then 

oo 

Sj(V) < J2 J p(t,X,x)dt. 

Theorem 3.3. // t~' r p(t,x,x)dt < oo for some 7 > ; then 



00 



t J p(t,x,x)dt. 



x&X 

V(x) 



(Note that here, the process x(t) may not be transient). 

The following two results are valid in both transient and recurrent cases. These results 
are based on the method of partial annihilation, proposed in [3],[1]. In the discrete 
situation it is equivalent to the rank-one perturbation technique. 
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Consider, for a fixed xq G X, the process x(t) with the condition of annihilation at Xq. 
The corresponding transition probability pi(t,x,y) is given by 

dp i 

— = H jn, x,y^x , Pi(t,x ,y) = 0\ pi(0,x,y) = 6 y {x). (11) 
As easy to see, / °° pi(t, x, x)dt < oo. 

Theorem 3.4. (CLR estimate, the general case). For any a, a > and some C\{o~), 

oo 

N (V) < l + #{x : V(x) > a} + Cl (a) ^ V ^ J Pi{t,x,x)dt. 



Theorem 3.5. [LT estimates, the general case). The following two estimates hold for 
each a > and some c(a) > 0: 

1 r°° 
SJV) < A 7 + ^ -V V l+ \x) / Pl (t,x,x)dt, (12) 

SJV) < + ^^y"V(x) f r' f p 1 (t,x,x)dt. (13) 

Here A is the largest eigenvalue of H . 

Note that Theorem 13.41 not only covers the recurrent case, but also provides a better 
result than Theorems I3.1[ 13.21 in the transient case when the operator H = H a depends 
on a parameter a which approaches a threshold a = ocq, where the process becomes 
recurrent. In Theorems 13.11 13T2"} the integrals in t blow up when a approaches «o whereas 
they remain bounded in Theorem 13.41 A similar remark is valid for Theorem 13.51 where 
the threshold depends on the values of a and 7. 

In the case where a = 0, [lj contains a more detailed description of the results obtained 
in Theorems 13.1113.51 



4 The perturbation of the hierarchical Laplacian. 

Theorems I3.1H3.3I and Proposition 12. 8[ when applied to the operator (TTUT) . lead to the 
same bounds on N (V) and S y (V) as in the case of the standard Schrodinger operator in 
]R d with the dimension d replaced by the spectral dimension s/,. An essential difference is 
that, while d must be an integer, the spectral dimension Sh can be an arbitrary positive 
number. The corresponding bounds hold if s > 2, where s = Sh in the estimate on Nq(V) 
and s = 7 + ^- in the estimates on 5 7 (V). The right-hand sides in these estimates blow 
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up when s j. 2 (the integrals in t diverge when s = 2). For example, Theorem 13.11 with 
a = and Proposition 12.81 imply a usual estimate: 

No(V) < #{x E X : V(x) > a} + E ^ k/2 (s)> 2 < s ft < A 

The case s < 2 is covered by Theorems 13.41 13.51 In fact, these theorems are valid for 
any s > and the estimates proven there are (locally) uniform in s. Hence they provide 
a better result in the transient case s > 2 than do Theorems 13. Itl3~3l when s 4- 2, see [1]. 

In order to apply Theorems 13.41 13~5| one needs to know an estimate on p\ as t — > oo and 
both the annihilation point Xo and x are arbitrary. If a = 0, then only the integral J °° p\dt 
is needed, not pi itself. The corresponding results can be found in [4] (we concentrated 
on N (V) in [3], but S y (V) can be studied similarly). Theorem 13.41 with o = implies jl] 
the following Bargmann type result: 

N (V)<l + #{x: V(x)>l} + C 1 (s h ) £ \/(a;)p(xo,a;) 2 -^, Sft < 2, (14) 

x:V(ie)<1 

with Ci(sft) — 7- oo as — >■ 2. A more accurate estimate of J °° pidt leads [4J to estimates 
on Nq(V) for all Sh and with a uniformly bounded constant: 

Theorem 4.1. If e < Sh < £ -1 , Sh ^ 2, then 

N (V)<l + #{x: V(x)>l} + C 2 (e) £ V(x) [1 ± f^^{ " 1 , (15) 

s:V(a;)<l ^ v 7 ^ 

// = 2, i/ien 

iV (V) < 1 + #{x : y(x)>l} + C 2 £ V(x) Hl+ /[ X °' X)] . 

s In -f 

x:V(z)<l VP 

In this section, we will obtain an estimate for pi as t — > oo, which allows one to use 
Theorems 13.41 13.51 with arbitrary o > 0. We will restrict ourselves to the case where 
Sh < 2 and provide an estimate only on Nq(V). The following refined Bargmann type 
estimate is an immediate consequence of Theorem 13.41 and Proposition 14.31 which will be 
proven below: 

Theorem 4.2. If Sh < 2, then 

N (V)<l + #{x: V{x)>l} + C 1 {s h ) V 2 -^{x)[l + p 2 {x^x)f- s K 

x:V(x)<\ 

We will conclude this section with a proof of the estimate on pi as t — > oo. This 
estimate is needed to justify the refined Bargmann estimate stated above and to prove 
similar estimates for SL. 
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Proposition 4.3. The following estimate is valid 

Pi(t, x, x) < C — — , t > 1, p = p{x , x), a = 1 - — . 

Remark. We expect that, in the case of fractal lattices similar to the Sierpincki 
lattice, the same estimate will be valid for a random walk with annihilation at a point. 
Proof. Consider the function 

POO 

R x 1) (x,y) = / e~ Xt Pl {t,x,y)dt. (16) 
Jo 

It is well defined when Re A > and understood in the sense of analytic continuation for 
complex A G C + = {A G C : |argA| < 3n/4}. From ffTTl) it follows that satisfies 

(Ah - X)R\\x,y) = -5 y (x), x,y ^ x , R x 1} (x ,y) = 0. 

Hence R^\x,y) = R\(x,y) + cR x (x } xq) } which together with the second relation in the 
formula above implies that 

R{ x\ x ,y) = R\{x,y) ~ ^ X Xo,V R x (x,x ). 

rC\{x , Xq) 

We put here y = x and -Ra^o? x) — R\( x o, x ) + -Ra^O; x) where (see Proposition 12.51) 

Rx(x ,x) = - r ~^~~^^2 f\ i s^ 8 ' r = d h (x Q ,x). (17) 

(\+p r l )v T v ^(A+p s )z/ s 

Taking also into account that R\(x,Xq) = R x (xo,x) and R x (x,x) does not depend on x, 
we obtain that _ 



R\ (x,x) = —2R\(xq } x) 



R x (x ,x )' 

We note that (1171) immediately implies the following two estimates: 



\R\{x ,x)\ < J^yi \R\{x ,x) - Rq(x ,x)\ < for all A G C+, r > 0, 

which together with (1181) and the Remark after Proposition 12.81 lead to 

R[ 1 \x,x) = a(r)+g(\,r), a(r) = -2R x (x , x), \g\ < + (19) 



(pis) r (pi/) 



The last estimate is valid for all A G C + with |A| < 1 and all r > 0. 
Applying the inverse Laplace transform to ( fTB"l) we obtain 



i rb+ioo 

Pi(t,x,x) = — e xt R i x \x,x)d\, 6»1. 

2vr Jfe.joo 
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Since is analytic in A G C + , and \R^\ < |j^r (the resolvent does not exceed the 
inverse distance from the spectrum), the last integral can be rewritten as 

Px(t,x,x) — — e xt R^\x,x)d\, 

where T = dC + with the direction on T such that ImA increases along T. We now use 
(|T9|) . the decay of R^ on T at infinity, and the fact that j r e xt d\ = 0, t > 0. This leads 
to 

/ x 1 f , / 2c|A| Ci|A| a \,,,, ai a 2 
Pl (t,x,x)<- / \e Xt \ 7 ^ + 7 ^)\d\\ = —±- + 



2tt ir \(P z/ ) r (P u ) 2r J t 2 (pv) r t 1 + a (pu) 2r ' 

It remains to recall that a = 1 — (see Proposition 12. 8j) . Thus pu = p a , and 



i — i — ^2 , 1 \Q 



(p 2 + l) a . □ 
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